A microscopic model of electronic field noise heating in ion traps 
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Motional heating of ions in micro-fabricated traps is one of the open challenges hindering experi- 
mental realizations of large-scale quantum processing devices. Recently a series of measurements of 
the heating rates in surface-electrode ion traps characterized their frequency, distance, and tempera- 
ture dependencies, but our understanding of the microscopic origin of this noise remains incomplete. 
In this work we develop a theoretical model for the electric field noise which is associated with a 
random distribution of adsorbed atoms on the trap electrode surface. By using first principle calcu- 
lations of the fluctuating dipole moments of the adsorbed atoms we evaluate the distance, frequency 
and temperature dependence of the resulting electric field fluctuation spectrum. Our theory repro- 
duces correctly the d~ 4 dependence with distance of the ion from the electrode surface and calculates 
the noise spectrum beyond the standard scenario of two-level fluctuators by incorporating all the 
relevant vibrational states. Our model predicts a regime of 1/f noise which commences at roughly 
the frequency of the fundamental phonon transition rate and a thermally activated noise spectrum 
which for higher temperatures exhibits a crossover as a function of frequency. 

PACS numbers: 37.10.Ty, 34.35. +a, 37.10. Rs, 72.70.+m 



INTRODUCTION 

Laser cooled trapped ions represent one of the most 
promising systems for the implementation of large scale 
quantum information processing jTJ [2] . Most of the basic 
requirements for building a quantum computer - the so- 
called DiVicenzo criteria [3] - have been demonstrated 
in the lab and the generation of entangled states of up 
to 14 ions has been achieved [I]. Many experimental ef- 
forts are now focused on the development of miniaturiza- 
tion and micro-fabrication techniques for ion traps [51-11 1 j , 
to realize more efficient and also fully scalable quantum 
computing architectures [Hill]- However, when devices 
are miniaturized, physics at the short-distance becomes 
a challenge. This is evident in measurements of Casimir 
force (of either sign) [13j [14] , or of non-contact friction 
[T51 - H7] , and in the case of trapped ions manifests itself in 
the appearance of an excess ("anomalous") heating rate 
as the trap-surface distance is decreased [Tl fTtJlHTlfTBH^] . 
Therefore, a detailed understanding of the origin of this 
noise will be essential for the future progress of trapped 
ion quantum computing, as well as the development of 
several hybrid quantum computing approaches where e.g. 
ions [221 IM| j Rydberg atoms [55] , polar molecules [25] or 
charged nano-mechanical resonators |27j are operated in 
the vicinity of solid state systems. 

Theoretical attempts to explain the noise-induced 
heating of trapped ions - its distance, frequency and tem- 
perature dependence - have been largely phenomcnolog- 
ical. The most common noise source in conductors, the 
Johnson-Nyquist noise from the trap electrodes or cir- 
cuitry, has a frequency independent spectrum and decays 
as dr 2 with increasing trap-surface distance d. Experi- 



mental data, however, is consistent with a d~~ 4 scaling for 
a large variety of trapping geometries (see e.g. Ref. [TO] 
for a recent review) and suggests a 1/f or even stronger 
variation over the observed range of frequencies [3 [20] . 
Therefore, since the early work of Turchette et al. [18j . 
the influence of randomly oriented dipole domains (patch 
potentials) has been recognized as the leading source for 
motional heating of ions. The main assumption of the 
patch potential model, namely that the electric noise 
originates from uncorrelated sources, explains correctly 
the observed distance dependence of the heating rate, 
but the model does otherwise not provide further insight 
into the physical origin of these fluctuations. More recent 
experiments with superconducting ion traps [19] strongly 
support the fact that the source of anomalous heating is 
not in the bulk, but is a surface phenomenon, and de- 
tailed temperature studies [5T] reveal that thermally 
activated processes are at play. 

In this work, we develop a microscopic theory of 
anomalous heating in ion traps and other charged sys- 
tems. The central assumption in this work is that the 
electric field noise in ion traps is produced by randomly- 
distributed fluctuating dipoles on the gold electrode sur- 
face, which in turn are formed by surface adsorption of 
atomic impurities, from the atmosphere or in fabrication. 
This assumption is in agreement with the experimental 
observations mentioned above and is supported by the 
fact that in many traps the heating rate increases over 
time, especially in the trap loading zone j 10] . 

In our model, fluctuations of the adatom dipole mo- 
ment arise from phonon induced transitions between mul- 
tiple bound surface states. We use analytic models sup- 
ported by exact density functional calculations to analyze 
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FIG. 1. (Color online) a) Fluctuating dipoles associated with 
a random distribution of adsorbed atoms lead to heating of an 
ion trapped at distance d above the surface, b) Dependence 
of the electric field fluctuation spectrum Se{uh) (in arbitrary 
units) as a function of the trap-surface distance d. The results 
have been obtained from a numerical simulation of the electric 
field from N = 100 uncorrelated dipoles distributed randomly 
over an area of 100 x 100 x do- do is the minimum separation 
between two dipoles. 



the adatom surface interaction potentials and the result- 
ing time variations of the induced dipole moment. From 
this analysis we obtain the electric field fluctuation spec- 
trum, and thereby the ion heating rate, as a function of 
the relevant microscopic parameters of the atom-surface 
interaction. Our calculations go beyond the standard sce- 
nario of two level fluctuators [7] 10, 28, 29] and we show 
that the inherent multi-level structure of the surface po- 
tential leads to a characteristic frequency and tempera- 
ture dependence, where even for a single atom a region 
with 1/f scaling emerges due to a distribution of differ- 
ent transition rates. The predicted distance dependence 
and heating rates are in good agreement with experimen- 
tal measurements. More importantly, the characteristic 
features of our noise model could provide more insight 
into the microscopic origin of anomalous heating and be 
tested against other potential mechanisms [30] . 



ANOMALOUS HEATING OF TRAPPED IONS 

Fig. [T] (a) shows a typical experimental setting where 
a single ion of mass m; and charge q is trapped at a 
distance d above a metal surface using combined static 
and inhomogeneous rf electric fields (Paul traps). For 



various designs of surface patterned micro traps [SHU], 
the separation d is between 20 to 100/xm. When the ion 
is laser cooled the vibrational ground state fluctuating 
electric fields couple to the motion of the ion and lead to 
an increase of the average vibrational occupation number 
n with a characteristic rate IT51 
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Here uit is the trapping frequency of the ion, typically 
in the range ui t /2n ps 0.1 — 10 MHz, and Se{u) = 
dT{5E{T)5E(0))e tu)T is the spectrum of the fluctuat- 
ing electric field SE(t) (projected onto the trapping axis) 
at the position of the trap. The heating rate ([lj can be 
measured in experiments and thereby provides an accu- 
rate probe of the electric field noise over the accessible 
distance and frequency range. 

In our model we consider electric fields originating from 
a distribution of fluctuating dipoles pi (t) which are asso- 
ciated with individual atoms adsorbed on the surface at 
positions fj. In accordance with previous "patch poten- 
tial" models [IS] [21] the assumption of uncorrelated noise 
sources leads to the expected distance dependence d~ 4 . 
This is illustrated in Fig. [I] (b) where we have numeri- 
cally evaluated the electric field noise of N = 100 ran- 
domly distributed dipoles on a surface. More explicitly, 
by averaging over a homogenous distribution of atoms 
with area density a the electric field noise spectrum at 
the position of the ion can be written as 



S>(w t ) 



8 (47re ) 2 d 4 



(2) 



where S^ui) — ^ oa dT(5^ z {T)5pL z {G))e lu)T is the spec- 
trum of the fluctuating dipole of a single adatom. Our 
main goal in the remainder of this paper is to provide a 
microscopic derivation of the dipole-fluctuation spectrum 
Sfj,(uj), which by using Eq. |l]) and Eq. pi) allows us to es- 
tablish a direct relation between the ion heating rate and 
the microscopic details of the atom surface interactions. 



ATOM-SURFACE INTERACTION 

Atoms approaching a surface experience an attractive 
force which at large distances is the well-studied van der 
Waals (vdW) potential ~ — C3/2 3 which eventually be- 
comes repulsive again when the electronic wavefunctions 
of the adsorbant and the bulk atoms overlap. An atom 
approaching the surface can lose energy by phonon in- 
duced processes and get trapped in the resulting potential 
well. We develop a model, based on physical intuition, 
which captures the essential aspects of the atom-surface 
interaction. To this end, we present ab initio density- 
functional calculations (DFT) of hydrogen adsorption on 
Au(lll) surface and obtain the interaction potential nor- 
mal to the surface. This potential is used as proxy for our 
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FIG. 2. (Color online) Potential energy curve for a H atom 
interacting with a 2 x 2 Au(lll) monolayer along the nor- 
mal to the surface calculated using density functional theory 
within the local density approximation. 
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FIG. 3. (Color online) Surface potential with potential depth 
Uo and minimum at zo- The potential support several bound 
vibrational states with energies hvi. The arrows represent the 
phonon mediated transitions between the bound states. The 
parameters used for this plot are Uo = 12 meV, zo = 6.05ao, 
P = 0.95a,7 1 and an atomic mass of m = 20 a.u. These 
parameters correspond to the case of Ne on a gold surface. 



model atom-surface potentials which more closely mimic 
realistic impurity reactivity on gold surfaces. The utility 
in using a parametric model potential rests in its flexibil- 
ity for tuning phonon transition frequencies, and range 
of both short-range and vdW interactions. 



Heavier adatoms on the gold surface mass scale the in- 
teraction potentials, leading to lowering of the transition 
frequencies. 



Model atom-surface potentials 



Ab initio atom-surface potentials 

For the calculations to be manageable, we chose hy- 
drogen adsorption on the gold surface. Clearly, surface 
contaminants are more than hydrogen atoms and contain 
atmospheric species. We incorporate such species by de- 
vising realistic interaction potentials [] to model the ad- 
sorption of atmospheric or fabrication impurities on gold 
surfaces. For the H-Au interaction, all-electron scalar 
relativistic calculations of the total energy and optimized 
geometries of a H-Au surface model system were per- 
formed using the spin-polarized density functional the- 
ory as implemented in the DMol3 software j32J. The ex- 
change correlation energy was calculated using the local 
gradient approximation (LDA) with the parametrization 
of Perdew and Wang (PWC) [33]. 

In Fig. [2j the H-Au(lll) surface potential is shown. 
The long-range interaction of atom-conductor potential is 
given by C3/2: 3 , where C3 is proportional to the atomic 
electric dipole transitions, and z is the normal to the 
surface. The coefficient C3 = J °° a(iuj)doj can be ob- 
tained numerically by evaluating the dynamic atomic po- 
larizability at imaginary frequencies. The value of C3 for 
hydrogen is given in Ref. [31] as C3 = 7.36 x 10~ 5 a.u. Al- 
lowance is made for the fact that Au dielectric constant 
is not infinite, i. e. Cz{Au) = fzjCa, where e is the 
dielectric constant for Au. The fundamental transition 
frequency is more than 25 THz, many orders of magni- 
tude larger than any frequency scale in the experiments. 



To describe the atom-surface interactions for wider 
range of atoms we use a model potential which is com- 
monly referred to as exp-3 [35 . This potential is fre- 
quently used in surface science studies [35] and provides 
a suitable description for short and medium range dis- 
tances. Our potential in this range is described by 
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where zq is the equilibrium position, Uq is the depth of 
the potential and /3 is the reciprocal range of repulsion. 
The typical shape of this potential is plotted in Fig. [3] In 
the long range Eq. (J3| reproduces the correct ~ — C3/2: 3 

form where C3 = gf^z3 Uo m terms of the model param- 
eters, and by adjusting these three parameters, we can 
modify U(z) to fit realistic atom-surface potentials. The 
potential supports several bound vibrational states with 
energies Ei = hvi as indicated in Fig. [3] By using an 
harmonic expansion of U(z) we find that for an adatom 
of mass m, the typical vibrational excitation frequency is 
approximately given by 



I Up 3(fi 2 - 4/3) 
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where j3 = fizQ. From this result, we also find a rough es- 
timate of the number of strongly bound vibrational states 
Af b w U /(hv 10 ). 
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From the exact H-Au potential shown in Fig. [2] we 
deduce Uq ~ 2 eV, corresponding to a temperature of 
T = 1.6 x 10 4 K, z w 1.6l and /3 w 3.91A" 1 , as ex- 
ample for an adatom with a high reactivity with the Au 
surface. Because of the low mass of H, also the vibra- 
tional frequencies are in the range of i/iq/2tt ps 40 THz. 
In general we expect a similar reactivity for other al- 
kali atoms as can be seen for example for K-Ag where 
U a = 1.79eV and z 2 A 36J. However, due to the 
larger mass we obtain significantly lower vibrational fre- 
quencies vio/2n 4 THz. In contrast, for weakly inter- 
acting atoms, we find from the widely-studied noble gas- 
metal potentials (see for instance, Refs. [38l HQ] ) that 
potentials can be wider and much shallower. For exam- 
ple, for Ne on a gold surface Uq « 0.012 eV, z s» 3.1A 
and /3 « 1.86A -1 . Within this typical range of potential 
parameters and adatom mass m ~ 10 — 150 a.u., we ex- 
pect the relevant vibrational frequency scales to be in the 
viq/2it w 0.1 - 1 THz regime. For our model potential 
in Fig. [3j and using the reduced mass of Ne-Au, we find 
v w /2tt « 0.3 THz. 

In our model of atom-surface interactions we ignore 
the dependence of U(z) on x and y due to surface rough- 
ness. However, for a metal surface, where the electrons 
are smeared out, this variation should be weak. Also, 
the random motion of an otherwise fixed dipole along the 
xyplane would lead to a different, d~ 6 scaling of the heat- 
ing rate which is not supported by experimental data. 

FLUCTUATING DIPOLES 
Adatom dipoles 

Adatoms adsorbed on the surface exhibit a finite per- 
manent dipole moment perpendicular to the surface. 
This induced dipole moment can be understood from the 
distortion of the electronic wavefunctions. It is commonly 
argued that the dipoles form, when the impurity adatom 
valence electrons penetrate into the surface conduction 
bands, modifying the surface work function. Here, we 
calculate the magnitude of the electric dipole moments 
of the adatoms, by resorting to elementary electrostat- 
ics, using image charge techniques. The electrons and 
the ionic core interact with the surface electrons through 
their respective image charges, as depicted in Fig. |4j We 
begin by writing the potential seen by the orbital electron 
due to its image charge [37], 

where fi = (pi, Zi) is the position of the i th electron and 
Z n is the distance from the nucleus of the adatom to the 
metallic surface (see Fig. Eh. In specifying the position 
of the electron, we use z, to refer to the distance from the 



electron to the nucleus along the normal to the surface. 
The first term is the above-discussed vdW interaction 
with the surface. The second term in Eq. ([5| vanishes in 
the first order of perturbation expansion, but contributes 
to the second order energy shift, 
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FIG. 4. (Color online) The adatom nucleus is at a distance 
Z n from the surface. The electron is at a distance fi from the 
nucleus. 

Following Ref. [37], we use a variational ansatz of 
the form <f> = (1 + AiZj)^>oj where Vo are the un- 
perturbed atomic wave functions, and calculate AE( 2 \ 
by varying the parameters A.;. From the deformed elec- 
tronic wavefunction <fi, we find an approximated induced 
dipole moment of the atom near the surface, P(Z n ) — 
(<j)\ J2i e f\4>)- For hydrogen atom P(Z n ) = 4 ' 5 ^ , where 
ao is the Bohr radius. Since the numerical factor in the 
expression is just the static polarizability of hydrogen, 
an = 4.5 aj], this expression can be generalized to a 
generic atom with polarizability a, 

P(Z n ) =0A7eal /2 a 3 / 2 ^. (6) 

We use the numerically constructed wave functions \i), 
corresponding to the vibrational bound states of our 
model potential U{z), to evaluate the average induced 
dipole moments 

^ := (i|P(*)|i). (7) 

We should note that due to the image charges the dipole 
moment seen by the atom will be twice the induced dipole 
moment in Eq. ([6]) . In Fig. [5] we plot the resulting in- 
duced dipole moments fj, Zi i for the model potential pa- 
rameters shown in Fig. [3] These parameters represent a 
weakly bound adatom similar to the Ne-Au surface po- 
tential [38]. Using this potential together with the po- 
larizability of Ne, a(Ne) = 0.361 3 [39 , and Z n = z , 
we obtain an induced ground state dipole moment of 
Mz,o(Ne) = 0.005D. Hence the dipole moment seen by 
the Ne atom is approximately 0.01 D which is in good 
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FIG. 5. (Color online) Magnitude of the average induced 
dipole moments p, Z: i for different vibrational states The 
values are shown for the model potential parameters given in 
Fig. [|] and a polarizability of a(Ne) w 0.36A 3 . 



agreement with 0.016 D calculated for Ne in 



Note 

that fi Zt i ~ a 3 / 2 /zg and in general the typical magni- 
tude of induced dipole moments is ~ ID. For example 
in [XT] HI] , it is shown that the dipole moment for Cs 
absorbed on Cu(100) is 4D. Similarly in [42] the induced 
dipole moments for K adsorbed on W, Ni and Pd ranges 
from 1.45 to 3.1 D. 



Phonon induced transition rates 

The emission and absorption of phonons in the bulk 
lead to transitions between different bound vibrational 
states and result in fluctuations of the induced dipole 
moment of the adatom. To evaluate the phonon induced 
transitions rates, we approximate the trap electrode sur- 
face by a semi-infinite crystalline slab. The position of 
an atom in the solid is given by 



u l5 



where r° is the equilibrium position of the i th atom and 
Ui is its deviation from the equilibrium. We can write 
the displacements terms of bosonic operators a\(q) 
for the phonon eigenmodes, 
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where q is the quasi-momentum, TV is the number of 
atoms in the bulk and M is their mass. For each q the 
normalized vectors e\(q) describe the three orthogonal 
phonon polarizations. 

In the presence of phonons, the adatom surface po- 
tential U(z) which in Eq. (3) has been defined for a 
static surface, will in general depend on the fluctuating 
positions of the bulk atoms so that U(z) — > U(z, {r^}). 



Expanding this potential to first order in gives, 

U(z, {r,}) ~ U(z, { ri }) + V VU(z, {r,°})u<. (9) 



Since our potential is already averaged over the two 
transverse directions, we are only interested in the vari- 
ations in the normal direction. Further, the dominant 
deformation of the potential arises from the closest sur- 
face atom and in the following, we restrict the sum in 
Eq. (U to a single atom. We use Fermi's Golden Rule to 
evaluate the phonon induced transition rate between two 
vibrational states \i) and |/), 
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where AE — Ei — Ej is the difference between the vi- 
brational energies Eij, n q are phonon mode occupation 
numbers and in this equation q = (q, A) includes the po- 
larization label. Assuming AE > and using the mode 
decomposition in Eq. (J8|, the resulting phonon emission 
and absorption rates can be written as 

^^\{f\j-U{z)\i)\M^f), (12) 



where Aw. 
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v is the averaged speed of sound in the surface material 
and p is its bulk density. 

In the following, we denote by To = Ti^(T = 0), 
the zero temperature decay rate from the first excited 
to the lowest vibrational state. From a simple harmonic 
approximation of U{z) around its minimum at z (il we 
obtain the scaling 
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Using v = 3962 m/s and p = 19.3 g/cm 3 for Au, and 
the model potential parameters given in Fig. (|]), with 
^io/2tt = 0.3 THz, we find r /27r 3.31 MHz. For K-Au 
where j/io/27r = 4 THz ro/27r is approximately 67 MHz 
and in general we expect r /27r to range from about 1 to 



a few hundred MHz. Note that the validity of Eq. (11) 



and Eq. (12) is restricted to transition frequencies Awj/ 
smaller than the Debye frequency ujd of the bulk mate- 
rial, which for gold is about 3.6 THz. 

Dipole fluctuation spectrum 

We are interested in the fluctuation spectrum of the 
induced dipole moment of a single adatom, defined as 



dr « Mz (tK(0)) - < M ,(0)) 
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Summarizing the results from the previous sections we 
can write the dipole moment operator as \x z = J2i Pz,iPu 
where pi = is the projection operator on the vi- 

brational level Therefore, for given values of Pi the 
dipole fluctuation spectrum can be related to the set of 
two-point correlation functions (pi(t)pj(t + r)} of the vi- 
brational populations. The populations in turn evolve 
according to the master equation 



dt 



(Pi) =^2M ij {p j 



(15) 



where the diagonal Ma = — Yljjti ^i^j an d the off- 
diagonal elements Mij = Tj^i are determined by the 
phonon induced transition rates discussed above. We 
evaluate the dynamics of the correlations (pi(t)pj(t + r)) 
by first introducing the condition Y^Api) = 1 the 
master equation. For N bound states we have, 



dt 



(Pi) = Al M + M ^(! - E 

= {M*j - M iN )( Pj ) + M iN , (16) 

Since we are only interested in level populations, all co- 
herences, pij = in the density matrix can be omit- 
ted. Using the quantum regression theorem, we find for 
* < N, 



d_ 
dt 



( Pi (T)p k (0)) = ^(M ij -M iJV )< ft (r)p fe (0))+M i ^ o) 

(17) 

where p k °^ is the steady state population in level k. For 
i = N we obtain 
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All two-point correlation functions can be calculated 
from the above two equations, and the full noise spectrum 
is obtained by summing all the two-point correlations. 



RESULTS AND DISCUSSION 

In Fig. |6j we plot the typical behavior of the dipole 
fluctuation spectrum S^u) as a function of frequency. 
The spectrum has been evaluated for the model potential 
parameters given in Fig. [3] and for different temperatures 
T. The lowest temperature fcsT = 0.02[/ corresponds 
to a situation where the thermal energy is smaller than 
the vibrational energy hisi ~ O.IUq and only transitions 
between the two lowest vibrational states contribute to 
dipole fluctuations. In this case, the dipole fluctuation 
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FIG. 6. (Color online) The dipole fluctuation spectrum 5 M (oi) 
is plotted for six different temperatures T (given in units of 
the vibrational excitation frequency vio) and the model po- 
tential parameters given in Fig. [3] and Fig [5] The frequency 
is scaled by the zero temperature transition rate To. The 
curves in black (fcsT /hvio = 0.2,0.3,0.4) correspond to tem- 
peratures where only the two lowest vibrational states are 
populated. At these temperatures the spectrum resembles 



that of a two level fluctuator (see Eq. ( 19 1 ) and its size in- 
creases with increasing temperature. For higher temperatures 
(shown in red, green and blue) more and more bound states 
are thermally occupied and an intermediate regime emerges 
where the dipole fluctuation spectrum exhibits 1// noise scal- 
ing. For all temperatures, the crossover from the flat (white 
noise) regime to a l// 2 noise or 1/f noise scaling occurs at 
around uj c w Yo(n(vio) + 1). 



spectrum resembles that of a thermally activated two 
level system [29], where at low temperatures 
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For frequencies below the fundamental phonon transi- 
tion rate, uj < Tq, the spectrum is frequency indepen- 
dent (white noise) while above r the scaling changes to 
1 /uj 2 as one would expect for a two-level system j^Hj • The 
noise is thermally activated with a characteristic temper- 
ature scale T = /u^o/fes, which depending on the atomic 
species ranges from few Kelvin to above room tempera- 
ture for light and tightly bound adatoms like hydrogen. 

When the temperature exceeds the characteristic vi- 
brational energy T > hi>io/kB, more and more vibra- 
tional levels are populated and contribute to the dipole 
fluctuations (curves in red, green and blue in Fig.[6|. An 
intermediate frequency regime appears where the noise 
spectrum exhibits a 1// scaling. We emphasis that in 
contrast to standard models for 1/f noise, based on 
a random set of two level fluctuators with varying pa- 
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rameters [251 US] , this l/f scaling occurs in our model 
even for a single dipole and emerges from a distribu- 
tion of different vibrational transitions rates Tif which 
contribute to the dynamics. We find that the approxi- 
mate crossover between the white noise and the 1// noise 
regime occurs at u s» Tq x {n{vw) + 1). To check the 
validity of our model, we compare our spectrum to the 
measured values. Using Eq. [2] we can relate to 
Se(w). Since oc /^ 2 , we note that in the 1// noise 
region, KT 11 < S^uj) < 1CT 7 D 2 /Hz. The value for 
r ~ 1 — 10 MHz, and the l/f behavior sets in at 
10 - 100 T which corresponds to 10 - 100 MHz. Us- 
ing the coverage fraction 6 = 0.1 corresponding to about 
cr - 10 18 m~ 2 , we find 3.2 x 10~ 8 < ujS e (u) < 0.0032 
V 2 /m 2 at do = 10 /im. The experimentally measured 
values range between KT 7 - 0.001 V 2 /m 2 [ID]. 
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Temperature dependence 

The temperature dependence of the ion heating noise is 
more succinctly displayed in Fig. [7] as a function of scaled 
temperature in units of ^io- At low frequencies, i.e. in the 
white noise regime, the fluctuations are thermally acti- 
vated, showing a peak at ksT '/hvio ~ 1, while for higher 
temperature they fall off again as ~ 1/T. The suppres- 
sion of low frequency fluctuations with increasing tem- 
perature is in principle expected from a single two level 
fluctuator. In some models based on multiple two level 
systems with a distribution of activation energies noise 
increases with increasing temperature over the whole fre- 
quency range |29j . In our model this behavior is recovered 
for frequencies within the "l/f region" of the spectrum, 
oj ss 20 x Tq , where the temperature dependence matches 
some of the experimental findings described in [7] and fits 
an Arrhenius curve F(T) = S T e- T °l T with the param- 
eters, St = 10" 13 D 2 /Hz and T Q « 0.242 U = 56.9 K. 
Finally, for very large frequencies, uj 3> Tq, i.e. in the 
l/f 2 regime, noise scales linearly with T. 

SUMMARY AND CONCLUSIONS 

In summary, we have developed a microscopic model 
for electric field noise generated by fluctuating dipoles 
associated with adatoms on a metallic surface. We have 
shown that phonon induced transitions between differ- 
ent bound vibrational states cause fluctuations of the in- 
duced dipole moments and generate electric field noise 
which can contribute to the anomalous heating observed 
in surface ion traps. The analysis presented in this work 
has been largely based on analytic model potentials for 
the atom-surface interactions and induced dipole mo- 
ments, allowing us to characterize the resulting field fluc- 
tuation spectrum for a wide range of atomic or molecular 
species in terms of a small set of microscopic parameters. 



FIG. 7. (Color online) Temperature dependence of the 
dipole fluctuation spectrum vs. normalized temperature 
(fesT '/Hvio)- The three curves show the temperature depen- 
dence evaluated for frequencies in the "white noise regime" 
(u) — > 0), the "1// regime" (u = 20 x r ) and the high fre- 
quency, "l/f 2 regime" (cj = 100 x To). The other parameters 
are the same as in Fig. [6] 

While more accurate predictions will depend on the 
atomic species, our model explains correctly the magni- 
tude and the d~ 4 scaling of the observed electric field 
noise. In contrast to standard models for l/f noise, 
which assume a distribution of two level fluctuators, our 
analysis predicts that for adatoms a transition from a flat 
to a l/f regime in the noise spectrum should occur only 
at finite temperatures and at frequencies above a typical 
phonon transition rate uj > T . From our estimates, we 
find that this transition rate is in the range, but slightly 
above the typical ion trapping frequencies. This suggest 
that either heavier or more loosely bound adatoms are 
responsible for the noise, or additional mechanism like 
multi-phonon transitions or dipole-dipole interactions are 
at play leading to the emergence of even lower fluctuation 
rates. 

Our model predicts several distinct features for the 
electric field noise spectrum which appear at character- 
istics frequency and temperature scales of the adatom 
surface interactions. In future experiments the develop- 
ment of new trap designs |43j could allow for more tar- 
geted search for these predictions, e.g. by probing sam- 
ples where a clean metal surface is contaminated with 
a single, pre-specified atomic species. To access the fre- 
quency regime of fO to a few fOO MHz, similar heat- 
ing experiments could be done with nano-mechanical res- 
onators or carbon nanotubes, for which accurate optical 
and electrical readout schemes have been developed (see 
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e.g. Ref. [33] and references therein). 
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